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the grid as well as the Kyoto protocol’s conditions and the 
accompanied EU directives, (c) to analyse the cost and revenue 
structure of the optimised power sector and identify the part due 
to greenhouse gas allowances, and (d) to investigate the share of 
electricity production based on renewable energy sources in the 
optimised power sector and calculate the expected CO emissions 
as a function of time. Among the milestones of the study was the 
stochastic correlation of the variables as well as the calculation of 
the optimal usage of the power plants contributing to electricity 
production. In order to accomplish the above milestones and 
objectives, some computational techniques had to be implemen- 
ted, to overcome the difficulties arising from nonlinear interaction 
of the various stochastic variables. The model has been imple- 
mented for the case study of the Greek Power Sector. However, 
other case studies are possible, provided that historical data of fuel 
and electricity prices are available and adjusted to those particular 
cases. 

The study is structured in brief as follows: in Section 2, a 
literature review is presented. In Section 3 an overview of the 
model is given while in Sections 4 and 5, details of the 
mathematical and computational model are provided. Section 6 
includes the results of the model’s implementation. Finally, in 
Section 7 a critical analysis is presented with suggestions on 
strategic interventions and policy issues related with deregulated 
electricity production. 


2. Related studies 


In the field of power generation investments, even recently, 
traditional cash-flow methods of analysis were used to address 
power sectors’ structure [4]. However, stochastic analysis experi- 
ence increasing acceptance among the energy-related scientific 
community. The answer to the question “invest or wait” is dealt in 
recent studies as well as in the present study. One may 
discriminate two separate classes of models: (i) those optimising 
policy interventions in future generation mix thus suggesting 
state-originating power licensing procedures [2,5] and (ii) those 
who suggest investment options for private energy-related 
businesses [6-8]. The present study should likely be classified to 
the first (i) category of modelling. 

Robust models are developed in the above mentioned studies of 
category i (optimised portfolios) although an optimal usage of 
power plants might be also incorporated, i.e. not only finding the 
required annual capacity investments but also calculating their 
optimal capacity factor depending on the generation mix over 
time. Recently, the uncertainties arising from climate changes 
seem to gain increasing focus [9,10]. These studies address the 
major issues of emissions trading that emerge in parallel with 
severe climate changes. Despite of the excellent analyses, either 
the CO, trading costs-revenues are not calculated inside the 
optimisation’s objective function or there is no power-portfolio’s 
optimisation included. In Refs. [11-13] methodological, organisa- 
tional and financial issues addressing the transition phase to 
deregulated electricity markets can be seen, but an analysis based 
on the environmental issues as well as emissions’ trading issues, is 
not included. Concerning power-portfolio optimisation [5,14] 
integrated multi-objective models have been developed. An 
elaborate consideration of heat recovery in different power 
generation methods is included, thus addressing the combined 
heat and power (CHP) industry too. Cost killing policies in different 
geographical regions are also addressed. In Ref. [15] discrete 
probabilistic optimisation models have been analysed for the 
evolution of stochastic variables. The models mentioned above, 
were treated by either dynamic or stochastic programming 
techniques. The SQP (Sequential Quadratic Programming) ap- 
proach of the present study is a consequence of the built-in 


interaction of the various contributing stochastic variables which 
in turn implied nonlinear considerations in the optimisation 
algorithm. 


3. Methodology 
3.1. An overview of the model 


The electricity market is assumed to be deregulated, as far as 
the electricity prices are concerned. Nonetheless, the annual upper 
bounds of licence provisions from each competing technology are 
assumed to be controlled by the state, thus reflecting the current 
situation in the Greek Power Sector. A random-walk procedure is 
used in order to forecast the future electricity demand and forward 
prices. The same stands for fuel prices as well as for CO, allowance 
prices. Appropriate stochastic differential equations (SDEs) (either 
Geometric Brownian Motion (GBM—[16]) or Mean Reverting (MR— 
[17])) are numerically solved using an Euler-Marujama solver. The 
uncertainty of the algorithm is considered by utilising a Monte- 
Carlo iterative sub-routine embedding that solver. The next step of 
the study is to reconstitute the future, in an optimised form: the 
previously forecasted data, feed an optimisation computational 
procedure which reproduces the optimal future structure of the 
power sector, i.e. what capacities from each type should be ordered 
yearly and what should their cumulative energy production be for 
the next years. Optimality is determined when the maximum Net 
Present Value (NPV) of the system is achieved while simulta- 
neously meeting the total power demand. It is important to note 
that the “system” mentioned above includes the state, the grid, the 
society and every power producer contributing to electricity 
production. The state in turn includes the regularity authorities 
and the system operator which are presented in Appendix A. The 
optimisation procedure comprises of a nonlinear algorithm, based 
on the solution of the Karush-Kuhn-Tucker (KKT) equations which 
are approximated with Sequential Quadratic Programming (SQP). 
Investment costs are calculated using learning curves of expertise- 
cost relationship in a similar vein to [2]. Emission’s trading costs 
and related revenues are embedded in the optimisation procedure 
as contributors to the objective function. Finally, the Kyoto 
protocol requirements and the related EU directives are also 
considered as basic problem constraints. 


3.2. Stochastic simulation 


The simulation of stochastic variables is required for represent- 
ing their evolution over time. Several options are potentially 
available: the class of ARMA (Auto Regressive Moving Average) 
models is mainly useful for stationary time-series simulations and 
the ARIMA (Auto Regressive Integrated Moving Average) models 
which are used for non-stationary time-series [18]. The GARCH 
(generalised autoregressive conditional heteroskedasticity) type of 
models [19,20] on the other hand is mainly indicated for time- 
series with extremely high volatility levels (noise). 

In the present article, futures contracts (i.e. using forward 
prices) are rather simulated instead of spot prices thus represent- 
ing average prices over a time period which do not contain 
information about the very short term variations of electricity 
prices (as stated in Ref. [21]). By utilising a probability distribution 
fitting tool, it has been validated that the time-averaged forward 
prices as well as the other stochastic variables samples are 
normally distributed. 


Assumption 1. Since normally distributed samples are required 
for GBM or MR simulations, it may be assumed that the time- 
averaged stochastic variables follow a random-walk process: (i) 
the Ito process (GBM stochastic differential equation), when the 
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variables’ history fail to define a long run mean value (electricity 
demand and electricity and fuel prices) and (ii) the Ornstein- 
Uhlenbeck (MR stochastic differential equation) in the opposite 
case (CO, allowance prices). In the latter case, the long run mean 
represents market equilibrium. Therefore, in the present article 
forward prices are modelled with a GBM process. Both processes 
(GBM/MR) may be considered as particular cases of ARMA (or 
ARIMA) models when dealing with discrete time and they are 
mathematically formulated with Eqs. (2*) and (3*), respectively, 
presented in Appendix B. 


Concerning forward electricity prices, Keppo and Lu [6] also 
adopted the above assumption (GBM representation of forward 
price evolution). In addition, they assumed that jump-processes 
are absorbed by averaging spot prices over time and using the 
resulting forward prices in futures contracts. Barlow [22] 
developed GBM and MR derivatives for the evolution of forward 
prices and their correlation with spot prices, accounting the 
spike absorption effects. Audet et al. [23] Laurikka [7] Laurikka 
and Koljonen [8] and Kumbaroglu et al. [2], have also adopted 
various types of GBM or MR models for the projection of forward 
and spot electricity prices as well as for the evolution of other 
stochastic commodities like fuel prices and electricity demand 
loads. 

In the present work, the Brownian differentials of the stochastic 
variables were correlated based on past data in an attempt to 
capture endogenous dynamics. A further analysis on the endoge- 
nous interaction of electricity prices with other stochastic 
variables has not been attempted. It would rather raise challenges 
in terms of computing time and input fitting, considering that a 
multiple investors problem is dealt with many individual power 
plants belonging to 10 discrete power generation technologies. 
Moreover, sufficiently detailed data would be needed reflecting the 
dynamics of price variations in relationship with capacity 
instalments, start-ups or shut downs for all the participating 
technologies. However such detailed data were not available. 
Nonetheless, the definite impact of an individual power investor’s 
actions on electricity prices may be minimised when multiple 
investors are considered and also, provided that the following two 
assumptions are additionally made: 


Assumption 2. Market equilibrium is considered 


Assumption 3. The investors’ profit is a direct function of electric- 
ity price 


As documented in Ref. [21], competitors’ actions may affect an 
investor’s decisions in two ways: (i) the option value of waiting 
may be reduced as this may force other competitors to enter the 
market prior to the first investor thus capturing its potential share 
on the market and (ii) the upper limit of profit distributions may be 
reduced since new plant instalments usually contribute to price 
dropping. Which effect is stronger depends on the characteristics 
of price modelling and on the other investors’ strategies. For 
multiple investor models - in which Assumptions 1, 2 and 3 are 
valid - the two effects are cancelled out and therefore the impact 
on electricity prices is minimised while the investment strategy 
remains the same as if other investors were not represented in the 
model (as proved by Dixit and Pindyck [24]). Thus, in the present 
work, being a multiple investors model, the endogenous modelling 
of electricity prices is solely based on the Assumptions 1, 2 and 3 
and on their statistical correlation with the other participating 
stochastic variables. They are all simulated through a Monte-Carlo 
process [16,25] which is comprised of two steps: (a) producing 
multiple SDE solutions using an Euler-Marujama solver and (b) 
averaging the multiple solution paths in order to supply the 
optimisation algorithm with the required input data. 


3.3. Optimisation 


The general scope of SQP numerical algorithm is to approximate 
an objective function with an easier one. This sub-problem can then 
be resolved and used as the basis of an iterative process. A 
characteristic of a large class of older methods was the transforma- 
tion of the constrained problem to a basic unconstrained problem by 
using a penalty function for constraints that are near or beyond the 
constraint boundary. The constrained problem is then solved using a 
sequence of parameterised unconstrained optimisations, which 
converge to the constrained problem. These methods are now 
considered relatively inefficient and have been replaced by methods 
that have focused on the solution of the Karush-Kuhn-Tucker (KKT) 
equations. The KKT equations are necessary conditions for 
optimality for a constrained optimisation problem. If the problem 
is a so-called convex programming problem, meaning that it has 
convex objective and constraint functions, then the KKT equations 
are both necessary and sufficient for a global solution point. 
Concerning the numerical solution of KKT equations, SQP methods 
are mostly utilised in nonlinear programming. Schittkowski [26] has 
implemented and tested an SQP version that outperforms every 
other tested method in terms of efficiency, accuracy, and percentage 
of successful solutions, over a large number of test problems. Based 
on previous works [27-29], the method calculates the Lagrangian 
factors of the KKT equations. Each iteration includes an approxima- 
tion of the Hessian matrix of the Lagrangian function using a quasi- 
Newton updating method. This is then used to generate a QP sub- 
problem whose solution is used to form a search direction for a line 
search procedure. An overview of the SQP method is found in [30- 
33]. During the last years, great advances have been made 
concerning the realisation and implementation of the above 
mentioned works. More recently, together with the development 
of highly efficient processors (CPUs), appropriate software have 
been available which are able to reach globally optimum solutions 
while maintaining high standards of accuracy as well as computa- 
tional efficiency [34]. 


3.4. Contribution of the study 


Apart from the results produced, which focused on the Greek 
Power Sector, the methodological contribution of the study may be 
summarised to the following combination of key-points: 


(1) Optimisation of the power sector’s future structure using a 
nonlinear SQP approach able to improve - on average - the 
future state of the power sector. 

(2) Impact assessment of CO; trading to business plans’ efficiency 
and state licensing policy, within the context of Kyoto 
requirements and EU directives concerning CO, emissions 
reduction. 

(3) Focus given on the stability of the national grid. 

(4) The power plants’ usage factor is not predetermined but it may 
be resulted from the optimisation procedure, allowing long 
term planning of operational intensity. 

(5) The Brownian differential of the stochastic variables are 
correlated based on their historical data, thus allowing their 
endogenous modelling; correlated uncertainties aim to 
improve the performance of analysis in business plans. 


4. The mathematical model 


4.1. Design considerations 


The essential idea of the model is the so-called “policy 
intervention” in business and society [24]. Optimal policy 
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intervention is attempted through stochastic analysis: the option 
to remain on existing generation mix competes with the option to 
order additional capacities based on emerging power production 
technologies. An objective function of the power sector’s NPV over 
time is created and an optimal point is depicted through an SQP 
optimisation algorithm. This point indicates the investment entry 
timing which maximises the aggregate system NPV but also 
determines the upper bound share of capacities of each plant 
category (i), allowed to be ordered (in a specific time). In other 
words the state will be able to organise its power licensing 
schedule (annually ordered capacities) in an optimal manner. The 
impact on the electricity generation industry is straightforward: 
let us suppose that the variables fluctuating stochastically and 
forecasted, like prices of electricity and fuel (“revenues” and 
“costs’”), are actually realised in the future. Also, the computation- 
ally resolved optimal licensing programmes and/or state directives 
for electricity production are supposedly violated time-wise. This 
means that some of the plants belonging to a specific plant 


Yr I v+lead;+li fe; 
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theoretically foreseen [24]. The entire process is a Monte-Carlo 
process comprised of two stages: (a) approximating the underlying 
multivariate process using an Euler-Marujama numerical solver of 
the vector-valued stochastic differential equations (described in 
Appendix B) and (b) averaging the multiple - log-normally 
distributed - solutions; these outputs feed the optimisation 
objective function. 


4.3. NPV optimisation and numerical algorithm 


As mentioned before, the model focuses on the maximisation of 
the system’s NPV. For this reason an objective function describing 
the system’s NPV has been created. 


4.3.1. Objective function and variables 

The mathematical formulation below represents the system’s 
NPV, considering cumulative annual incomes and expenses 
amortised over operational life-time: 


p(1 F r) E Lizvôizväm iai - 8760 


v=1 į¿=1 Z=v+lead; 
Yr I v+leadi+li fei 


v=1 i=1 z=v+leadi 
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NPV(Xi=1 p „Xiz, Oi-1, ET Oiz1) = Max 
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v=1 i=1 Z=v+lead; 
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ILL 
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category (i) might operate during a time period characterised by 
suboptimal conditions due to the predetermined “revenues” and 
“costs” evolution path (i.e. a plant type (j) of category (i) which 
might still operate while its fuel costs are very high). Consequently, 
the aggregate system-wise realised profits (NPV) would be 
suboptimal compared to the case of a generic compliance with 
optimal investment timing, strategy and directives, therefore 
leaving fewer chances for profit for the individual players. 

The maximisation of the NPV is equivalent to minimisation of 
the costs (CO2 trading, fuel mix, fixed and variable costs) while the 
revenues remain unaffected, since they are derived externally from 
the Monte-Carlo simulation of the stochastic variables, thus giving 
an element of stochastic optimisation to the model. The energy 
demand is the critical constraint that prevents the model from an 
unrealistic infinitesimal minimisation of costs, thus characterising 
the model as a demand-driven study. 


4.2. Random-walk procedure 


The random-walk approximation method used in this study 
allowed the simulation of generalised stochastic processes, and 
provided flexible simulation architecture. The procedure sup- 
ported nonlinear relationships which are commonly found in SDE 
simulations. In the present study the Brownian motion differ- 
entials were correlated based on historical data in an attempt to 
simulate the corresponding stochastic variables endogenously. 
Multiple - normally distributed - paths were then simulated thus 
considering the uncertainty of the model. The Brownian driven 
samples resulted to log-normally distributed output projections as 


C fuel, , (1 ae Tye hs Otitis - 8760 


C fom; z(1 F r) ae) Cem 


(1) 
Lizy0 8760. 160% .(4 + pe) 
izv%izvImi.i effi -( gi r) . Pco, 
Lizv9izv4midi * 8760 - eMmix - (1 +r) 7” - poo, 


The uncertainty is introduced for the stochastic variables p,, 
Cues Pco, Which follow the averaged projection paths generated 
through the Monte-Carlo process (described in Section 4.2). The 
€Mmix is the emissions factor from the current conventional 
generating mix (in tn CO2/MWh). 


CONV f 
co. 

5 LizvOizvOmiMli - 8760- orf 
i=1 

CONV 

5 LizvOizvaāmiđli -8760 

i=1 

CONV 


` l Emco, 
_ i=l 
CONV 


5 LizvOizvaāmiđli - 8760 


i=1 


(1a) 


The need for an automated algorithm based on a universal 
time-wise optimisation of NPV and not on successive time 
sweeps, generated a triple sum which feeds the computational 
procedure with the necessary data that correspond to the entire 
examined time period. An innovation of the model concerns the 
terms: 
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effi 


(+e. Pco, 


i=1 z=v-+lead; 
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and 


RENV-+leadj+li fe; 
5 5 Lizvôizvaām;iai - 8760 - eMmix ` qd ar jr ` Pco, 


i=1 z=v+lead; 


which represent the costs and revenues, respectively, arising from 
the CO» emissions trading. The costs correspond to the expenses of 
obtaining the required emission allowances for conventional 
power plants, while the revenues correspond to the incomes from 
trading the emission allowances generated by using renewable 
energy sources. The trading incomes result from multiplication of 
the renewable energy generated with the emissions factor of the 
current conventional generating mix and the allowance price, i.e. 
assuming that the renewable energy generated replaces energy 
that would otherwise be produced by the current conventional 
generating mix of the country. 
The variables of the objective function are: 


(i) the orders (capacities) realised in the year (v): Xiv (MW) for 
investing on new power plants of a specific technology and/or 
fuel type (i) and, 

(ii) the usage factor 6;,,, of the different types of installed power 
plants (i) operating during year (z) and ordered during year (v). 


It is noted that 6;,,, is requested only for the power generating 
units fired by expensive fuels (e.g. lignite, natural gas), otherwise 
(e.g. solar PV, wind-farms) it is imposed 6;,, = 1. This means that 
the power plants with zero fuel cost should be used whenever their 
energy source is available, thus reflecting real world practice while 
slightly reducing the variables and the computational cost. 

The investment costs CI;,,, of a power plant type (i) in the year 
(v) depend on technical advances arising from long periods of 
cumulative experience on construction of such power production 
units. This fact is expressed by the learning curves methodology 
which can be mathematically formulated as follows: 

z 2000 loga [1-b] 
e e vi (2) 


2000 
v=1960 Civ 


The investment expenses of year 2000 are used as a reference 
value. 

Two scenarios were considered, concerning the electricity 
demand: 


- Inthe first scenario (forecasted prices and demand or briefly FPD) 
the future loads are forecasted using historical data (the relative 
numerical approach is described in Section 4.2). 

- In the second (elasticity related prices and demand or briefly 
ELPD) the demand is related to its price with the following 
elasticity formula (similar to Ref. [2]): 


d: =a- p7" (3) 


This formula expresses a condition realised by the assumption: 
“Electricity demand depends only on its price and does not have 
any antagonists or counterparts impacting its price. Also, its 
elasticity (E) may be considered as constant”. Therefore in both 
scenarios, the total demand and the forward prices are related in an 
attempt to capture market dynamics endogenously. 

Excessive produced energy considerations: It is noted that energy 
produced in excess to the demand is a condition occasionally 
occurred during the numerical iterations of the optimisation 
algorithm until it converges to an optimal solution. Any excess 
energy produced is assumed to have zero value, as it cannot be 
exploited by the system. On the other hand, the production costs of 
the non-served energy, contribute (being a part of the expenses) in 


the NPV objective function and therefore they are always 
calculated (see Eq. (1)) thus reflecting the current situation in 
the Greek Power Sector. As a result, the optimisation algorithm 
tends to eliminate any excess energy production, by iteratively 
attempting to reduce the system costs 


4.3.2. The constraints of the optimisation problem 


(1) Logical constraints 


Usage factor (6;,,) values fluctuate in the space: 0 < 6;,, < 1 
(izv = 0 means a non-operational power plant— 6;,, = 1 means a 
maximum possible operational time). Also the annual orders for 
new power plants Xi, should be positive numbers: X;, >0, Vi,v. 


(2) Natural resource availability 


Some technologies of power production are constrained by 
limited natural resources and/or fuel reserves which determine the 
upper limit of the electricity generation potential. The wind 
generators, the hydro plants and the geothermal plants are 
indicative examples of such types of electricity generation whose 
resources present an upper potential limit. This constraint can be 
represented mathematically as follows: 


z+lead;+li fe; 
Xiv < dyn; Vi (4) 


v=z-+lead; 


(3) Kyoto protocol and EU Directive 2001/77/EC constraints 


Since the emissions trading has been engaged in the model’s 
objective function, additional constraints should be entered 
concerning the compliance with the Kyoto protocol requirements 
and the EU Directive 2001/77/EC. These constraints are obligatory 
for every participating country. Greece is required to meet the 
following targets: 


(3a) the share of renewable energy sources should be at least 20,1% 
of the total electricity production by the year 2010. In 
mathematical formulation: 


v+lead;+li fe; REN 
SChizvizvGm iA; -8760> (20, 1%) 
z=v+lead; i=1 
v+lead;+li fe;REN+CONV 
( 5 LizvOizvaām illi " 8760) Vz>2010 (5) 


Z=V+lead; i=1 


(3b) the share of renewable energy sources should be at least 30% 
of the total electricity production by the year 2020. In 
mathematical formulation: 


v+lead;+li fe; REN 
> Liz Wizutm Ani -8760> (30%) 
Z=v+lead; i=1 
v+lead;+li fe~REN+CONV 
( LizvOizvGmidti- 8760) Vz>2020 (6) 


z=V-+lead; i=1 


(4) Social constraints 


Meeting the demand target is a basic condition ensuring social 
acceptance of power production at a national level. In the present 
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study the total annual electricity production is imposed to be at 
least as much as the predicted annual electricity demand. This 
constraint is split in two discrete branches: 


(4a) Total demand target (energy): As stated before, two discrete 
scenarios (FPD and ELPD) have been considered in the model 
regarding the projection of demand. In both scenarios, the 
electricity generation (in terms of energy) should exceed the 
projected demand thus being able to serve the yearly demand 
aggregations: 


v+lead;+li fe;[REN+CONV] 
XO Lizv9izvGmii- 8760>d,-(1+ MAD) Vz (7) 


Z=v-+lead; i=1 


The rationale for this constraint was based on the assumption 
that a possible generic failure of power supply (system black-out) 
would result to an excessive social cost. Therefore this case should 
be avoided in the expense of additional capacity orders causing a 
system NPV reduction. An additional safety factor has been 
implemented, namely the forecast error. In the FPD scenario, the 
MAD (mean absolute deviation) indicator has been used as a 
measure of success of the forecast for a predetermined time period 
used for its validation. However, - given the increased uncertainly 
of the lognormal forecasted demand arrays -, this MAD value has 
been scaled up in order to adjust to the entire time horizon of the 
model. In the ELPD scenario a value derived from the elasticity 
formula and the electricity price forecast MAD has been used. 


(4b) Peak-power demand target: The constraint implies that the 
aggregate installed capacities should be able to serve the 
projected peak-power demand. The peak-power demand - 
annual past observations - follows a linear low-slope trend 
and therefore it has been linearly projected to the future. The 
availability and usage factors (except those of the wind 
turbines and the solar PVs) have been considered to be equal 
to 1 meaning full load operating conditions with available 
energy potential for the peak-power demand period. A safety 
factor (generation reserve margin taking values in the range 
0,1-0,2) has been further multiplied to the peak-power 
demand in order to account for a possible non-availability of 
wind or solar energy or the possibility of some plants’ failure: 


v-+Head;+li fe;[REN+CONV] 
Lizyv=Pp,-(1 +rm) Vz (8) 


Z=V-+lead; i=1 


Two discrete assumptions have been formed for each one of the 
demand scenarios concerning the peak-power demand estima- 
tions: 


(i) FPD scenario: The peak-power demand - as observed in each 
past year - followed a linear trend and therefore has been 
linearly projected to the future. 

(ii) ELPD scenario: This scenario generally predisposed a moderate 
growth concerning the demand uncertainty. Therefore, the 
peak-power demand has been related to the electricity prices 
in the same manner with the total electricity demand forecast 
(using Eq. (3)—elasticity related prices and demand). 


This constraint might be able to model any possible case of the 
system’s reliability limits (either total energy or peak power 
depending on which constraint proves to be stricter). The 
generation reserve margin concept provides additional safety on 
the model and has been also adopted in the works [2,5]. 


(5) Stability of the national grid 


It is well known that some power generation technologies 
which depend strongly on natural resources cannot be the base- 
load workhorses of electricity production. Such technologies are: 
the wind turbines, the solar PVs, the geothermic plants and in 
general most of the types of Power Generators that rely on 
renewable resources. Despite of their short setup periods and 
cheap-or-no need for fuels, they often suffer from resource 
unavailability and unpredictability and can severely affect the 
national grid’s stability. In the present study, it has been assumed 
that the cumulative annual electricity energy production based on 
such technologies should be less than 50% of the total annual 
demand: 


v+leadi+li fe; REN 
X Lizv9izv4midyi + 8760 < 50%-d, Vz (9) 


z=v-+lead; i=1 


(6) The limits of state licensing procedure 


The upper bound of the annually ordered capacities from each 
one technology should be less than 600 MWa, or X;, < 600 Vi, v. 
This value was a conclusion mainly based on the past experience 
for Greece, reflecting the limits of its economy’s robustness in 
conjunction with power demand. 


5. The computational procedure 


The data structures have been modified in order to create arrays 
with compatible time discretisation. Initially, the past data were 
available in different time intervals. Electricity prices (Euro/ 
MWh,,;) and demand (MW) were available in an hourly basis; CO3 
prices were available in a daily basis, while fuel prices and central 
bank rates were available in a monthly basis. This was 
contradictory with the requirement of supplying the optimisation 
algorithm with annual inputs. Compatible data structures were 
mandatory in order for the optimisation code to be usable. 
Therefore, prior to the optimisation algorithm runs, the historical 
data were subject to the following modifications: 


(a) Initially, some fuel prices, electricity prices and demand were 
missing. Linear interpolation techniques have been applied to 
overcome this problem. 

(b) The granularity selected for creating the history of stochastic 
variables was based on the widest available data sets [35] 
which for the available data were one month-periods. This 
means that e.g. daily prices can be averaged in monthly prices, 
but monthly prices could not be transformed in daily prices 
unless a regression or interpolation technique is applied in the 
expense of inaccurate variations modelling. Therefore the 
historical data were averaged in monthly basis unless they 
were already provided in that granularity. Concerning a 
forward electricity contract with maturity T1 and delivery 
period [T1, T2], an approximation of monthly averages F(T) 
(forward prices) can be realised if daily (hourly) prices f(t) are 
initially available within daily (or hourly) intervals dt according 
to the formula: 


1 T2 


EO) = 7T Jn 


(tdt Te [T1,T2] (10) 

It is noted that while the electricity prices were actually 
averaged, this was not the case for the demand loads. Demand 
forecasts were derived by multiplying the demand observa- 
tions d(t) by their occurrence time interval dt, and aggregating 
the products within each month thus resulting to total monthly 
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Table 1 
List of technology related parameters (source: Greek Ministry of Development). 
Integrated Oil Natural Lignite Biomass Solar Large Medium Hydro Geothermal 
coal gas gasifier PV onshore and high pumped 
gasification STAG (NH) wind turbines head hydro storage 
Investment cost (10? €/MWa) of the year 2000 1376 1127 1062 1888 2674 6555 1245 2491 3736 1350 
Fixed cost of year 2000 (operational, 70,80 SAT 2753 48,51 262,2 26,88 23,60 24,91 49,82 34,09 
maintenance, insurance, etc.) (in 10? €/MWg,) 
Availability factor 0,75 0,75 0,75 0,82 0,75 0,90 0,90 0,85 0,92 0,70 
Capacity factor 0,80 0,80 0,65 0,75 0,80 0,15 0,25 0,34 0,40 0,90 
Learning rate 0,05 0,01 0,01 0,01 0,15 0,20 0,10 0,00 0,00 0,00 
Commissioning time 4 3 3 4 3 2 1 10 10 3 
Efficiency factor 0,45 0,38 0,56 0,36 0,29 1 1 1 1 1 
Fuel CO, emissions (tn CO2/MWh fuel) 0,354 0,278 0,202 0,363 0 0 0 0 0 0 
Operational life-time (years) 30 30 30 40 20 30 30 45 45 20 
energy demand estimations (D(T) in MWh): the current mean average of Greece, while the interest rate, which 
T2 includes the inflation and a technological risk, is assumed to be 5%. 
D(T) = d(t)dt Te[T1,T2] (11) Another set of inputs necessary for the model was the initial status 


T1 


(c) When in monthly form, the variables were mathematically 
correlated using correlation factors derived by their past 
monthly history (either averaged or original) in a similar vein 
to [36]. Actually, the Brownian differentials of fuel, electricity 
prices and electricity demand were correlated (mutually as 
well as with the corresponding Brownian differentials of the 
interest and inflation rates). Especially concerning the histori- 
cal CO; allowance prices, they have been correlated with the 
historical electricity prices as has been adopted in Refs. [7,8]. 
The statistical parameters required for the GBM (or MR) 
simulation (drift, standard deviation, correlation factors) were 
extracted by the monthly arrays of the historical data. 

(d) The GBM (or MR) simulation of the stochastic variables 
produced multiple forecasts through the Monte-Carlo algo- 
rithm described in Section 4.2. These projections were then 
averaged in a yearly basis. Monthly or even shorter granula- 
rities have been also adopted in the past (e.g. [21]). Yearly 
discretised studies have been also conducted by [2,5]. 


The computational procedure used for the problem’s numerical 
modelling can be summarised to the following steps: 


Step 1. Determination of problem’s inputs and technological 
details, such as: setup periods, lead time (construction), 
operational life-time, basic investment and operational 
costs, and load factors. 

Step 2. Simultaneous prediction of electricity selling price, elec- 
tricity demand, CO, allowance and fuel prices. Correlated 
Brownian differentials - based on appropriately structured 
past data - are used through a multivariate Monte-Carlo 
process in which an Euler—-Marujama solver is used to 
produce multiple solutions of either GBM or MR modelling. 

Step 3. Supplying the optimisation algorithm with the numerical 
forecasts generated in step 2 (after appropriate post- 
processing of data structures). 

Step 4. Optimisation of the future power sector structure using an 
iterative SQP method. Reconstitution of future electricity 
production by plant type (and consequently by fuel source). 


5.1. The model’s input data 


Three basic parameters characterise the entire model: time, 
inflation and interest rates. The time period of the study is 50 years 
but only the first 25 years might include valid results as will be 
explained later on. The inflation rate is assumed to be 4%, close to 


of the Greek Electricity Generation system as of the year 2001 [10]. 
The entire field of technology related inputs needed for the 
iterative computational procedure is presented in table form 
(Table 1), including information for the various parameters and 
variables. 


5.2. Prediction of electricity prices and demand 


Among the various possible combinations of input parameters 
and assumptions tested with the model’s code, only two scenarios 
were selected for presentation. These two case studies represent 
reasonable hypotheses for the future of deregulated electricity 
markets: 


5.2.1. Case study A—explicitly forecasted electricity prices and 
demand (FPD) 

In deregulated markets the competition for domination usually 
pushes goods’ prices down. Electricity - being strongly inelastic - 
could be influenced by prices’ fluctuation only in rare circum- 
stances. Its evolution pattern usually follows the actual social 
needs through time. For this reason, the first case study was based 
on the assumption of increasing demand drift (slope), as observed 
in the historical data (step 1) of demand acquired by the Hellenic 
Transmission System’s Operator [37]. The electricity prices 
(national grid’s marginal prices) were also based on historical 
data acquired by the same source. This data actually constitutes 
the payments realised by HTSO S.A. to the electricity producers. 
The historical data used for both variables was initially available in 
an hourly basis for years 2001-2008, but they were computation- 
ally modified as described in the introduction of chapter 5. A 
Monte-Carlo procedure (multiple solutions of a GBM model as 
described in Section 4.2) produced the annual estimation of total 
demand loads and forward electricity prices which can be seen in 
Figs. 1 and 2, respectively. 

The standard deviation for the historical data of peak electricity 
demand has been calculated equal to 0,037 while its mean drift 
was 0,018. 

The stochastic differential standard deviation (defined by 
dividing standard deviation with mean stochastic differential) 
for the historical electricity prices has been calculated equal to 
14,97. 


5.2.2. Case study B—elasticity related horizon of electricity prices and 
demand (ELPD) 

The national grid marginal prices were extracted by HTSO S.A. 
as in the first case study. The same GBM model and the same 
Monte-Carlo numerical procedure were used in order to produce 
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Fig. 1. The forecasted electricity demand for FPD. 


their future field. The nominal electricity prices obtained were then 
transformed to real prices (constant prices with reference to the 
year 2001), based on the argument which implies that during an 
inflation period the expenses rise with almost the same rate as the 
incomes do, so prices of goods should be referenced to a constant 
time, in order to reveal their true future influence to demand, 
through elasticity. The prices’ reference to a constant time could be 
solely carried out using the assumed inflation rate, unless an 
investment cash-flow had been considered which would require 
interest rate-based transformation. 

A constant elasticity equal to 0,4 was assumed for the entire 
studied time period. The elasticity equation (3) was then applied in 
order to project the future electricity demand. It is noted that in 
order to compute the elasticity formula’s base, a simple algebraic 
equation was solved that related the price and the demand of 
electricity (both known values) during the first year of examina- 
tion. The assumptions made above were based on the initial 
hypothesis of non-existence of electricity substitutes impacting its 
price. This method finally resulted to a moderate increase of total 
demand as can be seen in Fig. 3. The forecasted forward prices 
retain the same path as in the previous case (FPD). Average 
electricity prices have been used for deriving the total demand 
(including peak loads). 

In both cases the expected demand increases with time but its 
slope is steeper in the FPD case. In the ELPD case the demanded 
loads are moderately increasing as a result of the constant 
elasticity assumption. Both cases comply with a reasonable 
hypothesis: “electricity demand keeps up with economies’ 
robustness and financial growth” [4,38]. It is noted that Greece 
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Fig. 3. The forecasted electricity demand for ELPD. 


experienced for almost a decade a constant economic growth in the 
levels of 3% or higher [39]. 


5.3. CO2 allowance prices 


The CO, allowance prices contribute to the system’s expenses 
and revenues as seen in Eq. (1). These prices present a constant - 
but slow - reverting to a long run mean derived by historical data 
[40]. As a consequence, the forecasting of these prices was 
performed by simulating an (MR) type of SDE (Eq. (3*) of 
Appendix B). The statistical analysis indicated a medium correla- 
tion with electricity prices (correlation factor 0,62). The forecasted 
CO, allowance prices can be seen in Fig. 4. 


5.4. Fuel prices 


Fuel prices contribute significantly to the finances of a power 
plant. Therefore, they should be forecasted paying particular 
attention on the correlation of their evolution patterns. This was 
done by simulating the GBM type SDEs (Eq. (2*) of Appendix B). The 
solution was repeated for almost 500 trials within the frame of an 
iterative Monte-Carlo simulation using correlated stochastic 
differentials as explained in the description of the mathematical 
model. 

The historical data of fuel prices cover the time period 2001- 
2008 and provide the necessary information about the statistical 
parameters required by the GBM simulation (mean drift, standard 
deviation, correlation factors, etc.). Past data sources were among 
others: the Hellenic Ministry of Development (hard copies) as well 
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Fig. 2. The forecasted forward electricity prices. 
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Fig. 4. The forecasted CO, allowance prices. 


Table 2 


Stochastic differentials’ correlation coefficient of fuel prices extracted from historical values. 


Hardcoal Oil Natural gas Lignite Biomass Geothermal Hydroelectric 
Hardcoal 1 0,87 0,89 —0,11 —0,57 1,00 0,78 
Oil 0,87 1 0,67 —0,20 —0,30 0,88 0,90 
Natural Gas 0,89 0,67 1 —0,39 —0,88 0,89 0,48 
Lignite —0,11 —0,20 —0,39 1 0,57 —0,13 0,13 
Biomass —0,57 —0,30 —0,88 0,57 1 —0,57 —0,04 
Geothermal 1,00 0,88 0,89 —0,13 —0,57 1 0,78 
Hydroelectric 0,78 0,90 0,48 0,13 —0,04 0,78 1 


as Ref. [41]. Contractual lignite prices were available and acquired 
from private mines. Concerning the biomass price, the technique 
established in Ref. [42] was applied, which is a Reverse Logistics 
approach augmented by a Holistic Activity Based Costing proce- 
dure. Having as a starting point the price calculated with this 
technique for the year 2001, a mean drift equal to the inflation rate 
was further applied in order to produce the biomass price 
evolution till the year 2008 (last year of input data). The table 
of stochastic differential correlations derived by the previously 
analysed historical data of fuel prices is presented (Table 2). 

The resulting fuel price evolution is presented in the following 
graph (Fig. 5) as an outcome of the GBM SDE solution augmented 
by the Monte-Carlo simulation. It is emphasised that the fuel prices 
include their mining and/or harvesting costs, their logistical costs 
(transportation, handling and warehousing) as well as the costs for 
their production or processing to their final form. 

The prices presented in the above diagram (values given in €/ 
MWh) follow in general their historically observed trend. For 


example oil prices are slowly reverting to a constant equilibrium 
which is a typical behaviour for oil products [24]. The anticipated 
coal prices are decreasing, which is in line with studies concerning 
the global coal market [43]. On the other hand, anticipated lignite 
prices present an increasing trend. The gradually depleting Greek 
reserves, might actually lead to harder winning efforts, thus 
justifying this behaviour. The statistical analysis indicated a 
definite correlation of lignite and electricity prices (0,77). Gas 
and electricity prices presented a slightly weaker correlation factor 
0,71. The previously analysed and processed data are gathered as 
inputs for the optimisation algorithm. 


5.5. The time periods of the study 


Within the frame of this study one could easily establish three 
discrete time periods: 

The initial aim of the study was to calculate the optimum 
capacity orders for the next 25 years, meaning till the year 2033. 
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Fig. 5. The forecasted fuel prices. 
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Fig. 6. The investment under increasing uncertainties. 
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However, the iterative optimisation procedure should run for a 
wider timeframe because some technologies require significant 
construction times (up to 10 years). Therefore, power loads and 
their financial influence should be calculated at least till the year 
2043. Another set of 7 years were included in order to let the 
computational code eliminate the time-boundary effects. For this 
reason, the upper time limit was imposed to reach the year 2050. 
Despite of true operational life-times, used for each power 
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producing technology, in order to extract more accurate results 
one should modify the inputs by selecting a much wider time 
period (i.e. till the year 2080), thus being able to investigate the 
influence of power plants to be built in the distant future. However, 
adopting an even wider timeframe would lead to two main 
obstacles: (i) In its current form the computational code examines 
the 50-year period and its CPU runtime approaches the limits of 
the acceptable. (ii) A wider time frame would actually cancel the 
validity of the random-walk forecasting procedure. Since its nature 
implies a constantly increasing uncertainty of the log-normally 
distributed results [24], it can be easily understood why the time 
limits of implementation should be tight (Fig. 6). Therefore, 
although the code runs for 50 years, only the results concerning the 
time period 2009-2033 may be considered valid, as was initially 
planned. 


6. Results of the code 


The results of the model are presented in the following sub- 
sections. 


6.1. The FPD case 
The capacity additions appear to have a complex profile as a 


function of time. In the valid range of results there is no plant 
category expected to dominate compared to others. Almost all 
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Fig. 10. The expected CO emissions per produced MWh, for FPD. 


plant types should be ordered to contribute in the increased power 
demands assumed for the FPD case. The expected yearly capacity 
orders can be seen in Fig. 7. 

The expected usage factor smoothly oscillates in high levels for 
the same reason (Fig. 8). 

In Fig. 9 the cumulative yearly electricity production is shown. 
It is interesting to note that electricity production succeeds to meet 
the demand target shown in the upper edge of the curves with a 
line. 


An interesting decrease of CO2 emissions per MWhg, produced 
is expected and this can be seen in the graph of Fig. 10. The 
anticipated increasing penetration of renewable energy sources - 
shown in Fig. 11 - might justify this behaviour. This is also an 
indication that the optimisation procedure attempted to lower the 
emission costs thus improving the system’s economy (i.e. 
compared to the current situation of year 2009). 

The computational code succeeded to follow the required 
increase of renewable energy penetration against the conventional 
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Fig. 11. The expected penetration of renewable energy (including hydroelectricity) in the Greek Power Sector. 


15 


A.I. Tolis et al./Renewable and Sustainable Energy Reviews 14 (2010) 2529-2546 


—*— Fixed Costs (investment etc.) 
—— Variable Costs (fuels etc.) 
~~-CO2 expenses 


2541 


EUROS-CURRENT PRICES 
E 


a 


25 30 35 40 45 50 


Fig. 12. The calculated annual costs and their expected breakdown. 


power plants. More specifically, starting from the levels of 2001 
(close to 8%) the renewables’ contribution rises to the current 
levels (close to 15% in 2008), and also manages to comply with the 
EU directives (20,1% in 2010 and 30% in 2020). 

The calculated costs breakdown is presented in Fig. 12. The 
expected share of CO, expenses - due to emissions trading - is 
moderate and slightly rising with time. However its percentage in 
comparison with the expected fuel costs (almost half) is significant 
enough to impact the oncoming investment decisions. 

The expected financial balance is presented in Fig. 13 for the 
entire system. The NPV of the balance curve corresponds to 
approximately 5 x 10° € as can be seen in Fig. 20. 
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6.2. The ELPD case 


The ELPD case produced a different decomposition of the Greek 
electricity sector. Within the range of valid forecasts (years 2009- 
2033, Fig. 14) hydroelectric and large onshore wind turbines 
should dominate the new orders, to maintain the medium levels of 
electricity demand arising from the constant elasticity hypothesis. 
Solar PV plants should be installed in the near future whereas the 
main base-load technology will remain the existing lignite plants. 

As can be seen in Fig. 15 the usage factor should vary in lower 
levels comparing to those observed in the FPD case especially for 
the lignite fired plants. This happens mainly due to the increasing 
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Fig. 13. The expected system's financial balancing for FPD. 
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Fig. 14. The calculated capacity orders for ELPD. 


costs forecasted for the case of lignite and the large share of lignite 
power plants at the starting period. 

High 90% are optimal levels for usage factors for all but lignite 
fired plants which should operate at 80% or lower for most of the 
time. The slight electricity demand increase of ELPD case is shown 
in Fig. 16 in comparison with the cumulative electricity produced: 
an obviously easier target to meet compared to the steep electricity 
demand increase of the FPD case. 

Due to the fact that the required additions are minor, the 
penetration of renewable energy sources is slower than in the FPD 


case. However the expected projection of emissions remains in 
significantly low levels too. As a result the expected CO, emissions/ 
MWh, produced are decreasing, as shown in Fig. 17 thus 
indicating a consequent economy improvement due to the 
optimisation algorithm. 

The expected expenses and financial balancing breakdowns are 
shown in Figs. 18 and 19, respectively. The CO, trading costs still 
contribute significantly to the system costs. 

It has to be noted that in Figs. 18 and 19, current (nominal) 
prices are presented meaning that their present value is 
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Fig. 16. The cumulative energy produced (ELPD) in comparison with forecasted demand. 
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Fig. 17. The expected CO, emissions per produced MWhg, for ELPD. 


significantly lower, thus resulting to a final NPV close to 7 x 10° € 
as shown in Fig. 20. 


6.3. Sensitivity analysis 


Despite of the model’s robustness for most of the input 
parameters, the system’s aggregate NPV could be severely 
influenced by the assumption of the interest rates. Therefore, 
the model’s sensitivity should be further investigated for a wider 
range of rates assumed. The computational procedure has been 


successively performed for acceptable rates in the range 4-8% with 
increasing step equal to 1%. The resulting sensitivity is shown in 
Fig. 20 for both scenarios. 

Lower interest rate values lead to higher NPVs as reasonably 
expected while a more moderate decrease appears for interest 
rates exceeding 6%. However, negative aggregate NPV is noted for 
high interest rate values, alerting for deficit systemic operation. It 
is noted also that the ELPD scenario proves to be less sensitive 
compared to FPD mainly due to the significantly lower need for 
future investments and fuel cost. 
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Fig. 20. Sensitivity of the model for the ELPD and FPD cases. 


7. Conclusions 


A computational model has been created in order to calculate 
the optimal power generation capacity additions in the context of 
future uncertainties. The model incorporates the calculation of the 
optimal usage factor of the power plant as well as the impact of 
emissions trading to the finance of the system. The aggregate 
system NPV constitutes the performance criterion of the study. 
Focus has been shown on the electricity sector of Greece. Two 
discrete scenarios have been examined in detail. In the first one, 
the electricity demand was forecasted directly from historical data 
(FPD) while in the second, a predetermined elasticity defined the 
future path of the electricity demand (ELPD). 

In the FPD scenario, the anticipated generation mix consists of 
existing but proven technologies as well as new emerging 
technologies based on renewable energy sources. The steep 
increasing of the electricity demand imposes the requirement 
for power generation originating from all the available technolo- 
gies, but not uniformly distributed. The low-or-zero fixed and fuel 
cost of wind generators and hydroelectric plants leads to a 
significant percentage of the required capacities being based on 


these technologies. However, fossil-fueled power plants like lignite 
and natural gas fired plants, should retain their position as 
strategic base-load fuels, as is the current situation in Greece. The 
remaining renewable source technologies like solar PVs and 
biomass-fired plants are challenged to this race and they may not 
gain a higher part of power generation unless experience acquired 
in the future reduces their investment cost. The usage of most plant 
types should remain in high levels close to 90% in order to meet the 
demand target. The produced electricity - as calculated by the 
model -, could comply with the Kyoto protocol requirements, 
concerning their future penetration against conventional power 
plants: a share exceeding 20,1% and 30% of the total production for 
the years 2010 and 2020, respectively should be based on 
renewable energy sources. Additionally, the stability of the system 
is ensured due to the restriction of renewable energies far below 
the limit of 50% of electricity demand. 

In the ELPD scenario the demand target is moderate and 
therefore easier to be met by the optimisation algorithm. The low 
investment and fuel cost of hydroelectric plants favours their 
future exploitation. Additional plant types based on renewable 
sources of energy like Solar PVs and biomass-fired plants would be 
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needed in low percentages, since the demand increase is not steep. 
For the same reason, the usage factor of future and existing lignite 
plant types should be moderate and rather lower than those 
calculated in the FPD scenario. 

The model proved to be robust in most of the input 
parameters but sensitive on the imposed interest rates. Regard- 
ing the performance of the two investigated scenarios, the ELPD 
scenario appears to be more efficient in terms of the aggregate 
anticipated NPV. The moderate requirements for electricity 
demand (ELPD) allow a higher penetration of renewable energies 
as a need for compliance with the Kyoto directives. Moreover, 
their lower fuel and operational costs contribute to a higher 
performance in terms of system’s economy. On the other hand in 
the FPD scenario, higher electricity demand levels are expected 
thus imposing the requirement for intensive use and higher 
capacity orders from both renewable and conventional plant 
types. Given the restricted energetic potentials of the renewable 
sources, the conventional plants should cover higher demand 
percentages in the expense of the aggregate system’s economic 
performance, especially due to their higher investment and 
emission costs. In both case studies a significant improvement 
with time in the COz emissions per MWh produced is expected, 
compared to the currently recorded values, thus indicating an 
optimised pattern in CO, emissions and therefore to the 
corresponding emission trading costs. 

The results indicated that the model could become a useful tool 
for energy policy advisors and authorities, state regulators and 
decision makers. Potential private investors and electricity-related 
business planners could also benefit from the tool and its results 
since they may have an overview of the future of different 
investment options and an indication of the shift in financial 
attractiveness of the various technologies due to the learning curve 
effects. As long as they may engage into the energy market arena 
they could identify their optimal investing options. 


Appendix A. Current Greek legal status notes 


In Greece the Law 2773/99 constitutes the basic legal background. 
Two companies, the Regulatory Authority of Energy (RAE) as well as 
HTSO S.A. (or the so-called “Hellenic Transmission System Operator” 
or even simpler “Operator” ) have been created. HTSO is the company 
that handles the operation of the Hellenic Transmission System of 
Electric Energy. These two companies are the basic key players of the 
deregulated electricity market. The main key players of the Greek 
Electricity market are: 


- RAE is an independent authority that manages, suggests and 
promotes the existence of equal opportunities and fair competi- 
tion and awards the operation license to producers, providers 
and all the other players, related to the market. 

- HTSO S.A. is a company with a twofold role: The first role is the 
one being played by P.P.C. in respect to the Transmission 
System: it ensures the balancing between production and 
consumption and that the electric energy is provided in a 
reliable, safe and in terms of quality acceptable way in terms of 
quality. The second role of HTSO is to settle the market, in other 
words to act like an energy stock market that arranges on a 
daily basis who owns to whom. HTSO does not provide electric 
energy and all the basic exchanging relations that take place 
are bilateral ones between producers/providers and their 
customers. 

- The Public Power Company (P.P.C.) is one among a number of 
companies that will be created in the field of electric energy. A 
stock market pictorial description of the roles could be, P.P.C. is 


an admitted company, HTSO is the stock market and RAE is the 
Capital Market Committee. 


Appendix B. The mathematical background for the forecasting 
tool 


Let us consider the following general stochastic differential 
equation (SDE): 


dXt = F(t, Xt)dt + G(t,Xt)dWt (1*) 
Where: 


e X is a state vector of process variables (for example, fuel or 
electricity prices) to simulate. 

e Wis a Brownian motion vector whose differential oscillates in a 
range of values generated by normal distribution. 

e F is a vector-valued drift-rate function. 

e Gis a matrix-valued diffusion-rate function. 


The drift and diffusion rates, F and G, respectively, are general 
functions of a real-valued scalar sample time t and state vector Xt. The 
Eq. (1*) is useful in implementing derived classes that impose 
additional structure on the drift and diffusion-rate functions. The 
derived classes of models used in this study are the following. 

GBM: The Geometric Brownian Motion Class whose general 
equation is: 


Where: 


e X; is a state vector of process variables. 

e y(t) is the mean drift function derived from the historical data as 
a mean difference. 

e D(t, X+) is a diagonal matrix-valued function. Each diagonal 
element of D is the corresponding element of the state vector 
raised to the corresponding element of an exponent (&œ(t)), which 
is a vector-valued function: D(t,X2 Je 

e V(t) is a matrix-valued volatility rate function. 

e dW(t) is a Brownian motion vector (noise) differential which is 
equal to € vdt and € € N(0,1). 


MR: The Hull-White/Vasicek (HWV) Gaussian Diffusion—Mean 
Reverting Class: 


dX, = S(t)[L(t) — X;|dt + V(t)dW, (3*) 
Where: 


e X; is a state vector of process variables. 

e S(t) is the reverting speed function derived from the historical 
data. 

e L(t) is the expected long run mean function. 

e V(t) is a matrix-valued volatility rate function. 

e dW(t) is a Brownian motion vector (noise) differential which is 
equal to € vdt and ¢ € N(0,1). 
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